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Abstract 

We give the approximate analytic solutions of the Dirac equations for the Rosen- Morse potential 
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I. INTRODUCTION 



Within the framework of the Dirac equation the spin symmetry arises if the magnitude of 
the attractive scalar potential S{r) and repulsive vector potential are nearly equal, S{r) ~ 
V{r) in the nuclei {i.e., when the difference potential A(r) = V{r) — S{r) — Cg — constant). 
However, the pseudospin symmetry occurs if S{r) ~ —V{r) are nearly equal {i.e., when the 
sum potential S(r) = V{r) + S{r) = Cps = constant) [1-3]. The spin symmetry is relevant 
for mesons [4]. The pseudospin symmetry concept has been applied to many systems in 
nuclear physics and related areas [2-7] and used to explain features of deformed nuclei [8], 
the super-deformation [9] and to estabhsh an effective nuclear shell-model scheme [5,6,10]. 
The pseudospin symmetry introduced in nuclear theory refers to a quasi-degeneracy of the 
single-nucleon doublets and can be characterized with the non-relativistic quantum numbers 
(n, l,j = / + 1/2) and (n— 1, / + 2, j = /-I-3/2), where n, / and j are the single-nucleon radial, 
orbital and total angular momentum quantum numbers for a single particle, respectively 
[5,6]. The total angular momentum is given as j — I + s, where I — I + 1 is a pseudo- 
angular momentum and 's — 1/2 is a pseudospin angular momentum. In real nuclei, the 
pseudospin symmetry is only an approximation and the quality of approximation depends on 
the pseudo-centrifugal potential and pseudospin orbital potential [11]. Alhaidari et al. [12] 
investigated in detail the physical interpretation on the three-dimensional Dirac equation 
in the context of spin symmetry limitation A(r) = and pseudospin symmetry limitation 
E(r) = 0. 

Some authors have applied the spin and pseudospin symmetry on several physical po- 
tentials, such as the harmonic oscillator [12-15], the Woods-Saxon potential [16], the Morse 
potential [17,18], the Hulthen potential [19], the Eckart potential [20-22], the molecular 
diatomic three-parameter potential [23], the Poschl-Teller potential [24], the Rosen-Morse 
potential [25] and the generalized Morse potential [26]. 

The exact solutions of the Dirac equation for the exponential-type potentials are possible 
only for the s-wave {I — Q case). However, for Z-states an approximation scheme has to be 
used to deal with the centrifugal and pseudo-centrifugal terms. Many authors have used 
different methods to study the partially exactly solvable and exactly solvable Schrodinger, 
Klein-Gordon (KG) and Dirac equations in ID, 3-D and/or any D-dimensional cases for 
different potentials [27-39]. In the context of spatially-dependent mass, we have also used 
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and applied a recently proposed approximation scheme [40] for the centrifugal term to find a 
quasi-exact analytic bound-state solution of the radial KG equation with spatially-dependent 
effective mass for scalar and vector Hulthen potentials in any arbitrary dimension D and 
orbital angular momentum quantum number / within the framework of the NU method 
[40-42]. 

Another physical potential is the Rosen-Morse potential [43] expressed in the form 

V{r) — —Vi sec h^ar + V2 tanh ar, (1) 

where Vi and V2 denote the depth of the potential and a is the range of the potential. 
This potential is useful for describing interatomic interaction of the linear molecules and 
helpful for disscussing polyatomic vibration energies such as the vibration states of NH^ 
molecule [43] . It is shown that the Rosen-Morse potential and its PT-symmetric version are 
the special cases of the five-parameter exponential- type potential model [44,45]. The exact 
energy spectrum of the trigonometric Rosen-Morse potential has been investigated by using 
super symmetric and improved quantization rule methods [46,47]. 

Recently, many works have been done to solve the Dirac equation to obtain the energy 
equation and the two-component spinor wave functions. Jia et al. [48] employed an improved 
approximation scheme to deal with the pseudo-centrifugal term to solve the Dirac equation 
with the generalized Poschl-Teller potential for arbitrary spin-orbit quantum number k. 
Zhang et al. [49] solved the Dirac equation with equal Scarf-type scalar and vector potentials 
by the method of the supcrsymmetric quantum mechanics (SUSYQM), shape invariance 
approach and the alternative method. Zou et al. [50] solved the Dirac equation with 
equal Eckart scalar and vector potentials in terms of SUSYQM method, shape invariance 
approach and function analysis method. Wei and Dong [51] obtained approximately the 
analytical bound state solutions of the Dirac equation with the Manning-Rosen for arbitrary 
spin-orbit coupling quantum number k. Thylwe [52] presented the approach inspired by 
amplitude-phase method for analyzing the radial Dirac equation to calculate phase shifts by 
including the spin- and pseudo-spin symmetries of relativistic spectra. Alhaidari [53] solved 
Dirac equation by separation of variables in spherical coordinates for a large class of non- 
central electromagnetic potentials. Berkdemir and Sever [54] investigated systematically the 
pseudospin symmetry solution of the Dirac equation for spin 1/2 particles moving within 
the Kratzer potential connected with an angle-dependent potential. Alberto et al. [55] 
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concluded that the values of energy spectra may not depend on the spinor structure of the 
particle, i.e., whether one has a spin-1/2 or a spin-0 particle. Also, they showed that a 
spin-1/2 or a spin-0 particle with the same mass and subject to the same scalar S{r) and 
vector V{r) potentials of equal magnitude, i.e., S{r) = ±V{r), will have the same energy 
spectrum (isospectrality) , including both bound and scattering states. 

In the present paper, our aim is to study the analytic solutions of the Dirac equation 
for the Rosen- Morse potential with arbitrary spin-orbit quantum number k by using a new 
approximation to deal with the centrifugal term. However, we use the approximation given 
in Ref. [56] which is quite different from the ones used in our previous works [39,40,42], 

^ cx"^ d + ^e-Ir)^ where d = or d = j^. The approximation given in [56] is convenient 
for the Rosen-Morse type potential because one may propose a more reasonable physical 
wave functions for this system. Under the conditions of the spin symmetry S{r) ~ V{r) and 
pseudospin symmetry S{r) ~ —V{r), we investigate the bound state energy eigenvalues and 
corresponding upper and lower spinor wave functions in the framework of the NU method. 
We also show that the spin and pseudospin symmetry Dirac solutions can be reduced to the 
S{r) = V{r) and S{r) = —V{r) in the cases of exact spin symmetry limitation A(r) = 
and pseudospin symmetry limitation E(r) = 0, respectively. Furthermore, the solutions 
obtained for the Dirac equation can be easily reduced to the Schrodinger solutions when the 
appropriate map of parameters is used. 

The paper is structured as follows: In Sect. 2, we outhne the NU method. Section 3 
is devoted to the analytic bound state solutions of the (3 -|- 1) -dimensional Dirac equation 
for the Rosen-Morse quantum system obtained by means of the NU method. The spin 
symmetry and pseudospin symmetry solutions are investigated. In Sect. 4, we study the 
cases AC = ±1 (/ = / = 0, i.e., s-wave), the Eckart-type potential, the PT-symmetric Rosen- 
Morse potential. Finally, the relevant conclusions are given in Sect. 5. 

II. NU METHOD 

The NU method [41] is briefly outlined here. It was proposed to solve the second-order 
differential equation of hypergeometric-type: 

<ir) + ^^Ur) + = 0, (2) 
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where a{r) and 5(r) are polynomials, at most, of second-degree, and T(r) is a first-degree 
polynomial. In order to find a particular solution for Eq. (2), let us decompose the wave- 
function ipni''^) as follows: 

^nir) = Hr)yn{r), (3) 

and use 

[<j(r)p(r)]' = T(r)p(r), (4) 

to reduce Eq. (2) to the form 

^Wnir) + r{r)y'^{r) + A|/„(r) = 0, (5) 

with 

r(r) = r{r) + 27r(r), r'(r) < 0, (6) 

where the prime denotes the differentiation with respect to r. One is looking for a family of 
solutions corresponding to 

A = A„ = -nr'{r) - (n - 1) a"(r), n = 0, 1, 2, • • • . (7) 

The Unir) can be expressed in terms of the Rodrigues relation: 

B r/" 

where Bn is the normalization constant and the weight function p(r) is the solution of the 
differential equation (4). The other part of the wavefunction (3) must satisfy the following 
logarithmic equation 

<P'ir) ^ 7r(r) 

(f>{r) a{ry ^ ^ 

By defining 

k = X-7:'{r). (10) 

one obtains the polynomial 

<r) = \ {a\r) - ?(r)] ± ^1 [^'(r) - ?(r)f - ?(r) + MO, (H) 

where 7r(r) is a parameter at most of order 1. The expression under the square root sign in 
the above equation can be arranged as a polynomial of second order where its discriminant 
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is zero. Hence, an equation for k is being obtained. After solving such an equation, the k 
values are determined through the NU method. 

In this regard, we derive a parametric generalization version of the NU method valid for 
any solvable potential by the method. We begin by writting the hypergeometric equation 
in general parametric form as 

[r (C3 - c^r)f ^^(r) + [r (cg - c^r) (ci - Csr)] ^^(r) + (-^r^ + - Cs) ^n(^) = 0> (12) 
with 

f(r) = ci - C2r, (13) 

a{r) = r (c3 - C4r) , (14) 

a{r)^-C,r' + C2r-C3, (15) 

where the coefficients q {i = 1,2,3,4) and the analytic expressions (j = 1,2,3). Fur- 
thermore, in comparing Eq. (12) with the counterpart Eq. (2), one obtains the appropriate 
analytic polynomials, energy equation and wave functions together with the associated co- 
efficients expressed in general parameteric form as displayed in Appendix A. 

III. ANALYTIC SOLUTION OF THE DIRAC-RO SEN-MORSE PROBLEM 

In spherical coordinates, the Dirac equation for fermionic massive spin-| particles inter- 
acting with arbitrary scalar potential S{r) and the time-component V{r) of a four-vector 
potential can be expressed as [26,57-60] 

[ca • p + /3 (Mc^ + Sir)) + V(r) - E] V^„,(r) = 0, V^„,(r) = V^(r, 9, 0), (16) 

where E is the relativistic energy of the system, M is the mass of a particle, p = —iKV is 
the momentum operator, and a and (3 are 4x4 Dirac matrices, i.e., 

• (s -i) (i ;) • (° "oO ■ (o -i) ■ 

where I denotes the 2x2 identity matrix and CTj are the three-vector Pauli spin matrices. 
For a spherical symmetrical nuclei, the total angular momentum operator of the nuclei J and 
spin-orbit matrix operator K = —f3 (<t • L -|- I) commute with the Dirac Hamiltonian, where 
L is the orbital angular momentum operator. The spinor wavefunctions can be classified 
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according to the radial quantum number n and the spin-orbit quantum number n and can 
be written using the Pauh-Dirac representation in the following forms: 

where Fni^{r) and GnK{f) are the radial wave functions of the upper- and lower-spinor compo- 
nents, respectively and Yj^{9, 0) and Yj^(9, cf)) are the spherical harmonic functions coupled 
to the total angular momentum j and it's projection m on the z axis. The orbital and pseudo- 
orbital angular momentum quantum numbers for spin symmetry / and pseudospin symmetry 
/ refer to the upper- and lower-spinor components, respectively, for which + — k{k + 1) 
and l{l + 1) — k{k — 1). The quantum number k is related to the quantum numbers for 
spin symmetry I and pseudospin symmetry I as 

(/ + 1) = - (j + i) , (si/2,P3/2, eic), j = l + l, aligned spin (k < 0) , 
+/ = + {j + I) , (pi/2, d3/2, etc.), 3 = l-\, unaligned spin («; > 0) , 
and the quasi-degenerate doublet structure can be expressed in terms of a pseudospin angular 
momentum 's — 1/2 and pseudo-orbital angular momentum I which is defined as 

-1 = - {j + I) , (si/2, P3/2, etc.), j = I - 1/2, aligned spin {n < 0) , 
+ (T+ l) = + (i + i) , {d3/2, /5/2, etc.), j = T+ 1/2, unaligned spin {k > 0) , 

where k — ±1, ±2, • • • . For example, (lsi/2, 0^3/2) and (2p3/2, 1/5/2) can be considered as 
pseudospin doublets. 

Thus, the substitution of Eq. (18) into Eq. (16) leads to the following two radial coupled 
Dirac equations for the spinor components 

{t + 7) ^"'^^''^ " ^^""^ + ^"'^ " ^^''^^ ^^^^^ 

- 7) Gn.{r) = (Mc^ - + E(r)) F„.(r), (19b) 

where A(r) = V{r) — S{r) and E(r) = V{r) + S{r) are the difference and sum potentials, 
respectively. 

Under the spin symmetry ( i.e., A(r) — Cg — constant), one can eliminate GnK{f) in 
Eq. (19a), with the aid of Eq. (19b), to obtain a second-order differential equation for the 
upper-spinor component as follows [16,26]: 
(P k{k + 1) 1 
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^ {El,-M\^ + Cs{Mc^-E^^))F^^{r), (20) 



where k{k + 1) = I {I + 1) , k = I for k < and k = — (/ + 1) for n > 0. The spin symmetry 
energy eigenvalues depend on n and k, i.e., — E{n, k{k + 1)). For I ^ 0, the states with 
j — I ± 1/2 are degenerate. Further, the lower-spinor component can be obtained from Eq. 
(19a) as 

where E^^ ^ —Mc^, only real positive energy states exist when Cg = (exact spin symme- 
try). 

On the other hand, under the pseudospin symmetry ( i.e., E(r) = Cpg — constant), one 
can eliminate ^^^(r) in Eq. (19b), with the aid of Eq. (19a), to obtain a second-order 
differential equation for the lower-spinor component as follows [16,26]: 
(P k{k-1) 1 



+ - ^ i^c^ - En. + C,.) A(r) 



= ^ - - i^c' + Gn.{r), (22) 

and the upper-spinor component Fni^{r) is obtained from Eq. (19b) as 

where Enn 7^ Mc^, only real negative energy states exist when Cps — (exact pseudospin 
symmetry). From the above equations, the energy eigenvalues depend on the quantum 
numbers n and k, and also the pseudo-orbital angular quantum number I according to 
k{k — 1) = /(/ + 1), which implies that j = I ±1/2 arc degenerate for / 7^ 0. The quantum 
condition is obtained from the finiteness of the solution at infinity and at the origin point, 
i.e., F„,(0) = ^^.(O) = and F„,(oo) = ^^.(oo) = 0. 

At this stage, we take the vector and scalar potentials in the form of Rosen-Morse poten- 
tial model (see Eq. (1)). Equations (20) and (22) can be solved exactly for k = 0, — 1 and 
K = 0, 1, respectively, because of the spin-orbit centrifugal and pseudo-centrifugal terms. 
Therefore, to find approximate solution for the radial Dirac equation with the Rosen-Morse 
potential, we have to use an approximation for the spin-orbit centrifugal term. For values 
of K that are not large and vibrations of the small amplitude about the minimum, Lu [56] 
has introduced an approximation to the centrifugal term near the minimum point r = re as 



1 1 



e 



_ _ — exp(— 2Q;r) _ / — exp(— 2ci;r) \^ 
1 + exp(— 2Q;rj Vl + exp(— 2Q;r) J 



(24) 
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where 

= -2 (exp(2.r.) + 1) [3 f il^^^^Iz^) _ (3 + 2.r.) ( ^ + -M'^ar.) 

2 /I + exp(— 2Q;re)\ ^ / 4Q;re 



D, = (eM2ar,)^ir^^^) ^3 + 2.., - ^^-^f^ ] , (25) 
and higher order terms are neglected. 



A. Spin symmetry solution of the Rosen-Morse Problem 

We take the sum potential in Eq. (20) as the Rosen-Morse potential model, i.e., 

(l + exp(-2ar))' ' (1 + exp(-2ar)) ^ ^ 

The choice of E(r) = 2V{r) V{r) as mentioned in Ref. [12] enables one to reduce the 
resulting relativistic solutions into their non-relativistic limit under appropriate transforma- 
tions. . 

Using the approximation given by Eq. (24) and introducing a new parameter change 
z{r) = — exp(— 2Q;r), this allows us to decompose the spin-symmetric Dirac equation (20) 
into the Schrodinger-like equation in the spherical coordinates for the upper-spinor compo- 
nent F„«,(r), 



dz^ z{l-z)dz z^l-z) 



2 



F^,{z) = 0, F„,(0) = F,,(l) = 0, (27) 



with 



/52 = ^ (2^0 + ^i) + ^ (Mc' + En. - C,) {Mc^ - E^^ - 21^)}, 



(28b) 
(28c) 



where uj — k{k-\-1) . 

In order to solve Eq. (27) by means of the NU method, we should compare it with Eq. 
(2) to obtain the following particular values for the parameters: 

t{z) ^1-z, a{z) ^zil-z), a{z) = -^^z^ + ^^z - s^- (29) 



Comparing Eqs. (13)-(15) with Eq. (29), we can easily obtain the coefficients Cj {i — 
1, 2, 3, 4) and the analytic expressions {j — 1, 2, 3). However, the values of the coefficients 
Ci {i = 5, 6, • • • , 16) are found from the relations A1-A5 of Appendix A. Therefore, the 

specific values of the coefficients Cj {i = 1,2, ■■■ ,16) together with C,j {j = 1,2,3) are 
displayed in Table 1. From the relations A6 and A7 of Appendix A together with the 
coefficients in Table 1, the selected forms of 7r{z) and k take the following particular values 



T^iz) = Snn " (1 + ^nn + S) Z, 



(30) 
(31) 



respectively, where 



(32) 



for bound state solutions. According to the NU method, the relations A8 and A9 of Appendix 
A give 

t{z) = 1 + 2enn - (3 + 2enn + 25) z, 

T'(r) = -(3 + 2£„, + 25) <0, (33) 

with prime denotes the derivative with respect to z. In addition, the relation AlO of Appendix 
A gives the energy equation for the Rosen-Morse potential in the Dirac theory as 



I 2fc2 2 

+ a n c 



iJhd^ + ~ + ^^^^ 



C^(gl+g2) 



- (n + 5 + 1) 



(n + 5 + 1) 

Further, for the exact spin symmetry case, V{r) — S{r) or Cs — 0, we obtain 

(Mc^ + E„,) {Mc' - E^, + V2) = -"^n^^ 



(34) 



I 2fc2 2 

+ a h c 



,2 



n 2 



with 



+ 5 + 1) 
5 = 5(C,^0). 



n + 5 + 1 



(35) 



(36) 
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Let us now find the corresponding wave functions for this model. Referring to Table 1 and 
the relations All and A 12 of Appendix A, we find the functions: 



p(z) = z^'"- (1 - zf+^ 

0(r) = (1 - zf^^ . 
Hence, the relation A13 of Appendix A gives 



(37) 
(38) 



y„(z)=^„^-2^"^(l-z)-(^^+^) 



1)^ 

dz"" 



n+2e„ 



(1 _ ^)"+2'5+l] ^ p^2e„.,25+l)(^ -2z), ze [0, 1], 

(39) 

where the Jacobi polynomial P!f''^\x) is defined only for ^ > — 1, z/ > — 1, and for the 
argument x G [—1, +1] • By using Fnniz) = 4>{z)yn{z), we get the radial upper-spinor wave 



functions from the relation A 14 as 



= Mnn (exp(-2ar))""'^ (1 - exp(-2ar))^+' 
X (-n,n + 2(£„„ + (5 + l);2£„«; + l;-exp(-2Q;r)). 



(40) 



The above upper-spinor component satisfies the restriction condition for the bound states, 
i.e., (5 > and > 0. The normalization constants Nnn are calculated in Appendix B. 

Before presenting the corresponding lower-component Gnnij'), let us recall a recurrence 
relation of hypergeometric function, which is used to solve Eq. (21) and present the corre- 
sponding lower component GnK{r), 



d 
dz 



iFi {a;b;c;z) 



,Fi (a + l;6 + l;c+l;^). 



(41) 



with which the corresponding lower component ^^^(r) can be obtained as follows 



A/;„,(-exi)(-2n/-))^"" (l + exi)(-2rw))'^+i 



2n (() + 1) exp(-2«r) k 
^'^ (1 -I- exp(— 2Q;r)) r 



+^fn. 



X 2F1 (-n, n + 2 (£„^ + 6 + 1); 2£„„ + 1; - exp(-2ar)) 

2an [n + 2 (sn. + S + 1)] (- exp(-2Q;r))""-+' (1 + exp(-2Q;r))^+' 
(2£„„ + l)(Mc2 + K«-C,) 



X 2F1 ( -n + 1; n + 2 ( + + ^ ) ; 2 (£„« + 1) ; - exp(-2Q;r) ) , (42) 
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where Enn 7^ —Mc^ for exact spin symmetry. Here, it should be noted that the hyper ge- 
ometric series 2F1 {—n, n + 2 + + 1) ; 2£nK + 1; — exp{—2ar)) does not terminate for 
n = and thus does not diverge for all values of real parameters 5 and 

For Cs > Mc^ + Enn and E^^ < Mc^ + V2 or Cs < Mc^ + E^^ and E^^ > Mc^ + 1^2, we 
note that parameters given in Eq. (28a) turn to be imaginary, i.e., < in the s-state 
{k — —1). As a result, the condition of existing bound states are > and S > 0, that 
is to say, in the case of Cg < Mc^ + Enn and Enn < Mc^ + V2, bound-states do exist for 
some quantum number k such as the s-state {k — —1). Of course, if these conditions are 
satisfied for existing bound-states, the energy equation and wave functions are the same as 
these given in Eq. (34) and Eqs. (40) and (42). 

B. Pseudospin symmetry solution of the Rosen-Morse Problem 

Now taking the difference potential in Eq. (22) as the Rosen- Morse potential model, i.e., 

A(r) = -W, exp(-2ar) (1 - exp(-2c.r)) 

(l + exp(-2Q;r))' ' (1 + exp(-2Q;r)) ' ^ ^ 

leads us to obtain a Schrodinger-like equation for the lower-spinor component GnniT)-! 



dz'^^ z{l-z)dz^ z^{l-zy 
where 



Gr^^iz) = 0, ^^^(O) = Gr,.{l) = 0, (44) 



(45a) 

^1 = ^ I ^ (i^o + /^i + /^2) - ^ [El, - M^c^ - {Mc" + K«) Cps - {Mc^ - En. + Cps) V2] 

(45b) 

(45c) 

and u; — k{k — 1). To avoid repetition in the solution of Eq. (44), a first inspection 
for the relationship between the present set of parameters (SnK, Pi, P2) and the previous 
set {EriK, Pi, P2) tells us that the negative energy solution for pseudospin symmetry, where 
S{r) — —V{r), can be obtained directly from those of the positive energy solution above for 
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spin symmetry using the parameter map [57-59]: 

Fn.{r) ^ Gnn{r), V{r) ^ -V{r) (or V, ^ -V, and V2 ^ -V2), ^ -En. and ^ -C^^- 

(46) 

Following the previous results with the above transformations, we finally arrive at the energy 
equation 



where 



(Mc' - K« + Cps) {Mc^ + En. - V2) = -^/I'c' 



2fe2 2 

a n c 



{n + Si + 1) 



- (n + 5i + 1) 



(47) 



(48) 



By using Gnniz) — (f){z)yn{z), we get the radial lower-spinor wave functions as 

Gnnir) = Mnn (exp(-2ar)f '"^ (1 - exp(-2ar))'^^+^ p^2?„«,25i+i)(^^ - 2 exp(-2ar)). (49) 

The above upper-spinor component satisfies the restriction condition for the bound states, 
i.e., 5i > and > 0. The normalization constants Hni are calculated in Appendix B. 



IV. DISCUSSIONS 



In this section, wc arc going to study four special cases of the energy eigenvalues given by 
Eqs. (34) and (47) for the spin and pseudospin symmetry, respectively First, let us study 
s-wave case I — Q {k — —1) and I — Q {k — 1) case 



(Mc' + En,-i - Cs) {Mc^ - En,-i + V2) = a^ffc^ 



2a'/i-c- 



n + 62 + l 



n + 62+l 



(50) 



where 



(61) 



If one sets Cg — into Eq. (50) and Cpg — into Eq. (47), we obtain for spin and pseudospin 
symmetric Dirac theory, 

J^(Mc2 + K,-i) 



(Mc^ + £;„,_i) (Mc^ - En,-i + V2) = a^h^c^ 



2a2n2c2 



n + + 1 



+ n + 5^1 + 1 



(52) 
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+ JI^<*^=' + ^"-''l• 



(53) 



and 



(Mc^ - K,+i) (Mc^ + - V2) = +a2/i'c2 



n + 5+1 + 1 



+ 71 + 5+1 + 1 



54) 
(55) 



respectively. The above solutions for the s-wave are found to be identical for spin and 
pseudospin cases S{r) — V{r) and S{r) — —V{r), respectively. 

Second, when we set Vi —Vi and V2 — > —V2 the potential reduces to the Eckart-type 
potential and energy eigenvalues are given by 

2 



(Mc' + (Mc' - En,-i - V2) = a'n'c' 



n + 5_i + 1 



+ n + (5_i + 1 



(56) 
(57) 



for spin symmetry and 



n + 5+1 + 1 



+ 71 + 5+1 + 1 



^..4 -l + Vl + ^(Mc-E„,+i)|, 



(58) 
(59) 



for pseudospin symmetry. 

Third, let us now discuss the non-rclativistic limit of the energy eigenvalues and wave 
functions of our solution. If we take Cg — and put S{r) — V{r) — 5](r), the non-relativistic 
limit of energy equation (34) and wave functions (40) under the following appropriate trans- 
formations (Mc^ + E^^) /h^c^ 2ii/h^ and Mc^ - E^^ -E^i [26,57,42] become 



Eni = V2 + - — - —a 



2i^ri 2// 



^ (21^1 + V2)-^^ + {n + ly + (271 + 1) 60 

7l + 5o+ 1 



with 



(60) 



(61) 
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and the associated wave functions are 

Fni{r) = N-ni {eM-'^o^r)Y- (1 - exp(-2ar))^+^° pi^e^Mi) _ ^ ^^^^_2o,r)), (62) 



where 



2a \l r? 



(63a) 



which are identical with Ref. [25] in the solution of the Schrodinger equation. Finally, the 
Jacobi polynomials can be expressed in terms of the hypergeometric function as 

pM(^l - 2exp(-2ar)) = ^Fi (-n, 1 + + i/ + n; + 1; exp(-2ar)) , (64) 

77/. 

where z G [0, 1] which lie within or on the boundary of the interval [—1, 1] . 

Fourth, if we choose V2 iV2, the potential becomes the PT-symmetric Rosen- Morse 
potential, where P denotes parity operator and T denotes time reversal. For a potential 
V{r), making the transformation of r — r (or r ^ ^ — r) and i —i, if we have the relation 
V{—r) — V*{r)), the potential V{r) is said to be PT-symmetric [43]. In this case we obtain 
for spin-symmetric Dirac equation 



I 2fc2 2 

-I- a Ti c 

■ n + 5 + 1 

In the non-relativistic limit, it turns to become 



-1 2 



- (^n + 5 + 1) 



(65) 



Enl = iV2 + 



7. a 



{2Vi + 1V2) -^ + {n + lY + (2n + 1) 5o 



(66) 



n + 5o + 1 

where real Vi > 0, which is identical to the results of Ref. [25]. If one sets Z = in the above 
equation, the result is identical with that of Refs. [44,45]. 



V. CONCLUSIONS 



We have obtained analytically the energy spectra and corresponding wave functions of 
the Dirac equation for the Rosen- Morse potential under the conditions of the spin symmetry 
and pseudospin symmetry in the context of the Nikiforov-Uvarov method. For any spin-orbit 
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coupling centrifugal term we have found the explicit expressions for energy eigenvalues 
and associated wave functions in closed form. The most stringent interesting result is that 
the present spin and pseudospin symmetry cases can be easily reduced to the KG solution 
once S{r) = V{r) and S{r) = —V{r) {i.e., Cg = Cpg = 0) [55]. The resulting solutions 
of the wave functions are being expressed in terms of the generalized Jacobi polynomials. 
Obviously, the relativistic solution can be reduced to it's non-relativistic limit by the choice 
of appropriate mapping transformations. Also, in case when spin-orbit quantum number 
K, — 0, the problem reduces to the s-wave solution. The s-wave Rosen-Morse, the Eckart- 
type potential, the PT-symmetric Rosen-Morse potential. and the non-relativistic cases are 
briefly studied. 
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APPENDIX A: PARAMETRIC GENERALIZATION OF THE NU METHOD 

Our systematical derivation holds for any potential form. 

(i) The relevant coefficients Cj (i = 5, 6, • • • , 16) are given as follows: 

1 1 

C5 = - (C3 - Ci) , C6 = - (C2 - 2C4) , C7 = Cg + ,^1, (Al) 
Cs = 2C5C6 - Cg = C5 + ^3, Cio = C4 (C3C8 + C4C9) + C3C7, (A2) 

2 2 
cii = — \/c^, C12 = \/cI^> (A3) 

C3 C3C4 

Cl3 = — (C5 + Vcg) , Cu = (\/cio " C4C5 - C3C6) , (A4) 

C3 C3C4 
2 1 

ci5 = — \/cio, C16 = — - C4C5 - C3C6) . (A5) 
C3 C3 

(ii) The analytic results for the key polynomials: 

7r(r) = C5 + - — (c4^ + ^/c^ - C3C6) r, (A6) 
C3 

= (C3C8 + 2C4C9 + 2VC9C10) , (A7) 
C3 

2 

T(r) = C3 + 2^ (C3C4 + C4^/c^ + ^/c^) r, (A8) 

C3 

2 

r'(r) = (C3C4 + C4^/c^ + v^) < 0. (A9) 

C3 

(iii) The energy equation: 

C2n-(2n + 1) C6 + — (2n + 1) (\/cI^ + C4\/c^)+n (n - 1) C4+\ (csCg + 2C4C9 + 2^/c^) = 0. 
C3 C3 

(AlO) 

(iv) The wave functions: 

p(r) =r'^"(c3-C4r)^l^ (All) 
0(r) = r^i3(c3 _ c^^)ci4^ > 0, ci4 > 0, (A12) 

yn«(r) = Pi^"'^^^Hc3 - 2c4r), Cn > -1, C12 > -1, r e [(03 - l)/2c4, (1 + 03)7204] , (A13) 
i^n.ir) = <P{r)yn.{r) = NnV^'^c^ - C4ry''PjC'''''''\c3 - 2c4r), (AM) 
where P^"'^^ (03 — 2c4r) are the Jacobi polynomials and J\fn is a normalizing factor. 
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When C4 — 0, the Jacobi polynomial turn to be the generalized Laguerre polynomial and 
the constants relevant to this polynomial change are 

lim Pj^^^^^'''^\c3 - 2c4r) = L^"(ci5r), (A15) 

C4— >0 

lim (cs — c^rY'^^ = exp(— Cigr), (A16) 
^n.{r) =A4exp(-ci6r)L^"(ci5r), (A17) 
where L^"(ci5r) are the generalized Laguerre polynomials and Afn is a normalizing constant. 



APPENDIX B: NORMALIZATION OF THE RADIAL WAVE FUNCTION 



In order to find the normalization factor AfnK, we start by writting the normalization 
condition: 



^ [' Z^'--\l - Zf^^'' [pP--,25+l)(i _ 2z)fdz = 1. (Bl) 

2« Jo 

Unfortunately, there is no formula available to calculate this key integration. Neveretheless, 
we can find the explicit normalization constant Mni- For this purpose, it is not difficult to 

obtain the results of the above integral by using the following formulas [61-64] 

1 

;i _ sf-^ s^-i ^F, {a, P; 7; as) ds = sF^ {ly, + T, «) , (B2) 



/ 



and 2F1 {a,b;c;z) = rui'irih) ^'""rfn+p^^^'' Hence, the normalization constants for the 



upper-spinor component are 



p=0 



T{26 + 3)r(25„, + 1) ^ (-1)"^ (n + 2(1 + + 6))^ T{n + m) 



2ar{n) 



m=0 



ml (m + 2ennV-T (m + 2 + 5 + 1)) 



-1/2 



, (B3) 



with 



fuK = 3F2 y^^uK + m, -n, n + 2(1 + + (5); m + 2 + 5 + - J ; 1 + 2£„„; ij , (B4) 

where (x)^ = r{x + m)/r(x). Also, the normalization constants for the lower-spinor com- 
ponent are 

-1/2 



T{26i + 3)r(2?„, + 1) (-1)"' {n + 2(1 + + ^^1))^ T{n + m) 



2ar{n) 



m=0 



m\ (m 



+ 2ennV.T {m + 2{enn + Si + l))' 



(B5) 
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with 



QuK = 3-^2 I 2£n« + m, -n, n + 2(1 + + m + 2 ( + 5i + - ) ; 1 + 2£„«; 1 ) . (B6) 
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TABLE I: The specific values for the parametric constants necessary for calculating the energy 
eigenvalues and eigenfunctions of the spin symmetry Dirac wave equation. 



Constant Analytic value Constant Analytic value 



Cl 


1 


C2 


1 


C3 


1 


C4 


1 







Cfi 
'^D 


1 
2 






C8 


-/32 


eg 


^riK 


ClO 




Cll 




C12 = Cl5 


2(5 + 1 


Cl3 




Cl4 = C16 


<5 + l 


^1 


Pi 


C2 


^2 


^3 


£2 
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